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This  report  concerns  the  design  of  a seismic  event  detection  assoc 
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cally convert  seismic  signal  detections  from  a ulIw  rk  of  stations  stored 
centrally  as  a list  of  station  detection  bulletins  into  a corresponding  event 
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gorithms process  single  site  short-period  detections  as  well  as  array 
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ABSTRACT 


This  report  concerns  the  design  of  a seismic  event  detection 
association  processor.  This  processor  consists  of  algorithms  which  auto- 
matically convert  seismic  signal  detections  from  a network  of  stations  stor- 
ed centrally  as  a list  of  station  detection  bulletins  into  a corresponding  event 
list,  while  discarding  spurious  detections,  or  false  alarms.  These  algori- 
thms process  single  site  short-period  detections  as  well  as  array  detections. 
An  automatic  depth  estimate  is  also  performed.  Testing  via  simulation  in- 
dicates that  the  resulting  event  list  is  accurate  enough  so  that  relevant  wave- 
form data  may  be  retrieved  automatically.  It  is  recommended  that  the  de- 
tection association  processor  be  further  evaluated  on  both' simulated  and 
real  networks. 


Neither  the  Advanced  Research  Projects  Agency  nor  the  Air  Force 
Technical  Applications  Center  will  be  responsible  for  information  contained 
herein  which  has  been  supplied  by  other  organizations  or  contractors,  and 
this  document  is  subject  to  later  revision  as  may  be  necessary.  The  views 
and  conclusions  presented  are  those  of  the  authors  and  should  not  be  inter- 
preted as  necessarily  representing  the  official  policies,  either  expressed  or 
implied,  of  the  Advanced  Research  Projects  Agency,  the  Air  Force  Technical 
Applications  Center,  or  the  US  Government. 
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SECTION  l 

INTRODUCTION 

The  detection  association  processor,  DAP,  is  the  basic  esti- 
mation and  control  system  for  an  automatic  seismic  network.  The  DAP  con- 
verts a stream  of  automatic  station  detection  bulletins  into  an  event  list  with 
enough  accuracy  to  retrieve  relevant  waveform  data.  The  DAP  can  operate 
on  an  array  network,  a single  site  network,  or  a mixed  network.  This  re- 
port is  primarily  concerned  with  the  documentation  of  those  algorithms  com- 
prising the  DA  P. 

The  report  consists  of  five  sections.  Section  11  defines  in 
greater  detail  the  role  of  the  DAP  in  an  automatic  seismic  network.  Section 
III  discusses  the  DAP  operation  for  an  array,  single  site,  and  a mixed  net- 
work. In  Section  IV  the  performance  of  the  DAP  upon  four  days  of  simulated 
network  data  is  analyzed.  In  Section  V conclusions  are  drawn  and  recom- 
mendations mcle  relevant  to  installing  the  DAP  as  an  automatic  netw  >rk 
element. 
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SECTION  II 


REQUIREMENTS  FOR  A DETECTION  ASSOCIATOR 

The  detection  association  processor  carries  out  the  function 
of  transforming  signal  detections  obtained  independently  by  a network  of  seis- 
mic stations  to  seismic  event  detections.  The  location,  depth,  and  origin 
time  of  the  estirm  ted  seismic  event  focus  are  then  used  to  edit  waveform 
segments  associated  with  the  events  by  estimating  arrival  times  of  the  seis- 
mic phases,  particularly  the  P-phase. 

The  requirement  for  information  supplied  by  the  detection 
bulletins  from  single  sensor  seismic  stations  is  the  seismic  station  identifica- 
tion, a detection  ranking  parameter  to  rank  bulletins  in  the  order  of  probabil- 
ity of  detection,  the  arrival  time  of  the  signal,  and  an  estimation  of  the  arrival 
time  error. 

Detection  ranking  parameter  depends  on  a station's  noise  sta- 
tistics. For  example,  if  a detector  sequentially  measures  , then  the  bul- 
letins may  be  ranked  by  m ; another  possibility  is  the  signal -to-noise  ratio, 

b 

S/N,  of  the  detections;  and  another  would  be  the  Z- statistic.  The  latter  nor- 
malizes a detector  output,  such  as  power  in  decibels,  by  subtracting  a run- 
ning estimate  of  the  mean  noise  level.  That  difference  is  then  divided  by  a 
running  estimate  of  the  standard  deviation  of  the  noise  level.  The  purpose 
of  7.  is  to  normalize  time-varying  noise  levels  and  noise  standard  deviations 
to  an  approximately  fixed  unit  normal  statistical  distribution.  Larger  posi- 
tive deviations  of  7.  , the  presumed  unit  normal  statistic,  are  related  to  low- 
er probability  of  a false  alarm  detection.  The  probability  of  the  false  alarm 
can  be  determined  from  normal  statistic  tables.  In  this  report,  Z- statistic 
estimates  were  used  as  the  detection  ranking  parameter. 
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For  array  seismic  stations,  additional  information  must  be 
provided.  This  is  direction  and  velocity  r(  the  signal  as  indicated  by  the 
maximum  beam  output,  and  the  corresponding  estimates  of  the  standard  de- 
viation of  those  quantities. 

The  detection  bulletin  list,  which  is  continuously  updated  in  real- 
time, is  input  data  to  the  DAP.  From  comparisons  of  arrival  time  and  other 
information  obtained  from  detection  bulletins,  preliminary  event  location  and 
depth  estimates  are  generated  in  real-time.  In  reporting  such  locations,  a 
time  lag  is  required  for  propagating  the  event  around  the  world  and  for  DAP 
processing.  The  time  lag  is  estimated  to  be  between  j hour  and  1 hour.  A 
time  requirement  for  event  locations  within  one  or  two  hours  of  real-time  is 
considered  reasonable  for  a DAP. 

The  purpose  of  the  DAP  is  to  support  those  seismic  analysts 
responsible  for  collecting  and  processing  segmented  event  records  associated 
with  seismic  events.  Once  the  analyst  obtains  records  of  event  signals,  the 
analyst  then  measures  accurate  arrival  times  of  selected  data  and  other  param- 
eters to  accurately  determine  the  event  focus,  magnitude,  and  other  event 
source  parameters.  Thus  it  is  seen  that  the  DAP  is  a tool  which  supports 
rapid  data  collection  by  seismic  analysts.  The  association  function  of  the  DAP 
is  to  automatically  associate  crude  arrival  time  measurements  and  other  meas- 
urements of  an  automatic  detector  with  all  events  detected  by  at  least  four  sta- 
tions. Given  one-to-one  association  of  four  or  more  arrival  times  with  an 
event,  a list  of  approximate  event  locations  can  be  generated  by  the  DAP. 

The  sole  operational  function  of  this  preliminary  seismic  event  list  is  to  either 
automatically  or  interactively  retrieve  seismic  event  edit  records  for  use  by 
the  seismic  analyst  who  routinely  performs  much  more  meaningful  evaluation 
and  measurement  of  the  possible  event.  The  preliminary  event  location  need 
only  be  accurate  enough  to  enable  the  analyst  to  request  records  of  some 
specified  length  assumed  to  be  considerably  greater  than  15  seconds  as  a 
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representation  of  desired  cha  rt - p«*  riod  seismic  phases.  For  this  data  re- 
trieval purpose,  locations  need  only  be  ai  c urate  to  about  2 at  teleseismit 
o 

distance  and  1 at  regional  distance. 

The  DAP  is  designed  to  merge  detection  bulletins  generated  by 
a world-wide  network  on  tin-  order  of  25  array  or  single  - sensor  stations.  The 
process  must  keep  up,  within  reporting  time  requirements,  with  a highly  vari- 
able rate  of  detection  bulletin  generation  by  the  network.  This  is  due  to  event 
swarms,  false  alarms  due  to  local  events  or  noise,  and  falsi-  alarms  due  to 
th«-  coda  and  later  phases  associated  with  larger  events.  It  is  estimated  that 
for  a 25  station  network,  a maximum  of  about  2000  bulletins  per  day  with  at 
least  10  words  per  bulletin  need  to  be  processed.  To  handle  this  volume  of 
information,  it  is  expected  that  most  of  the  DA  P ope  rations  need  to  be  per- 
formed in  an  automatic  mode,  possibly  with  quality  checking  and  station  and 
network  performance  monitoring  done  interactively  by  an  analyst.  In  this  re- 
port, designs  of  automated  DAP's  are  evaluated  using  simulated  detection 
bulletins  generated  by  an  earth  seismicity  model  and  a network  of  simulated 
station  detection  processors  given  by  Shoup  and  Sax  (1974). 

Another  requirement  of  the  DAP  is  to  obtain  algorithms  which 
can  perform  the  association  of  signal  detections  to  event  locations  for  a net- 
work of  single  sensor  stations,  a network  of  array  stations,  or  a mixed  net- 
work of  both  types  of  stations. 


SMC  Tlt.N  ill 


A SSC )(  I AT  OH  MODUL  MS 

The  re  .irr  llirt't*  major  associator  modules.  These  .1  it  the 
Arr.iy  Network  Processor,  Single-Site  Network  l’roeessor,  and  a Mixed  Net- 
work Processor.  The  Array  Network  Processor  operates  on  detection  bulletins 
from  automatic  detectors  at  array  stations.  The  Single-Site  Network  Process- 
or operates  on  detection  bulletins  from  automatic  detectors  at  sites  providing 
only  .1  single  vertical  component  measurement  of  signals.  The  Mixed  Network 
Processor  operates  on  detections  from  a network  of  both  types  of  stations. 

The  sub-system  functions  performed  by  the  network  are  as  follows: 

• Real-Time  Bulletin  Processing 

Add  new  bulletins  to  the  current  sta^  k of  detection  bulletins 
in  the  detection  bulletin  buffer 

Maintain  ranking  of  the  detection  bulletins  by  /.-statistic 
Remove  aged  bulletins  from  workspace  after  a fixed  num- 
ber of  failures  to  associate' 

Remove  bulletins  from  workspace  which  are  associated 
with  previously  located  events. 

• Association  Processing 

Set  up  one  or  more  preliminary  trial  locations 
Find  detections  which  are  consistent  with  the  locations 
Refine  the  trial  locations  and  repeat  the  association  tests 
Iterate  the  association  tests  until  a specified  number  of 
detections  are  associated  e.g.  tour  detections,  and  a 
local. ible  event  is  declared. 


Perform  final  location  of  a declared  locatable  event  with 
arrival  times  from  associated  direction  bulletins 


Use  depth  grid  and  depth  constrained  locations  to  find 
best  estimate  of  the  depth. 


Figure  III-l  illustrates  the  flow  of  decisions  and  actions  per- 
formed by  the  detection  association  processor.  The  inputs  are  detection  bul- 
letins  from  seismic  stations;  either  arrays  serviced  by  an  array  module  or 
single-site  stations  serviced  by  a single-site  module.  The  detection  associa- 
tion processor  operates  on  a fully  loaded  buffer  of  the  most  recently  acquired 
detection  bulletins  sent  by  seismic  station  detection  processors.  As  space 
becomes  available  in  tile  bulletin  buffer,  the  detection  bulletin  buffer  is  fed 
new  bulletins  from  a queue.  This  space  becomes  available  in  several  ways. 

One  way  is  to  fail  to  associate  with  different  events  some  specified  number  of 
times,  e.  g.  , three  times.  This  aging  process  is  a critical  factor  in  managing 
the  workspace  of  the  detection  bulletin  buffer.  If  the  aging  process  is  not  care- 
fully designed,  bulletins  will  age  out  because  the  workspace  is  clogged  with  coda 
detections  from  preceding  larger  events.  This  is  definitely  a problem  with  the 
HAP  models  described  in  this  report  and  this  problem  needs  to  be  attended  if 
better  results  are  to  be  obtained.  Another  way  for  a detection  bulletin  to  be 
discarded  is  to  have  been  associated  with  an  event  already  located  and  regis- 
tered on  the  event  list  generated  as  an  output  of  the  system.  Upon  discarding 
a bulletin  in  one  of  these  ways,  a new'  bulletin  is  fed  in  serially  from  the  queue 
and  inserted  into  the  detection  bulletin  buffer.  The  insertion  is  rank  ordered 
by  the  magnitude  of  the  detector  output,  so  that  the  most  probable  event  de- 
tections are  further  up  toward  the  top  of  the  stack.  First,  all  of  the  array  de- 
tections are  tested  for  association  with  the  presumed  location  of  the  top  ranked 
detection.  If  the  association  tests  are  passed,  the  location  is  updated  by  the 
Kalman  location  algorithm.  Otherwise,  a counter  is  incremented  to  register 
the  number  of  times  the  detection  bulletin  tailed  to  associate. 
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After  completing  the  association  tests  of  all  of  the  array  stations, 
tests  are  performed  on  the  single-site  detection  bulletins.  This  is  done  by  cen- 
tering a grid  of  potential  locations  on  the  location  estimate  obtained  by  the  pre- 
ceding operations  of  the  array  module.  If  there  are  no  array  detection  bullet- 
ins, then  the  location  of  the  station  with  the  largest  detector  value  is  taken  as 
an  initial  location,  about  which  a grid  of  potential  locations  is  centered.  See 
Figure  I II - 2 . Association  is  performed  by  the  single- site  module  by  search- 
ing for  the  maximum  number  of  origin  times  consistent  with  a depth  location 
at  each  grid  point  corresponding  to  a trial  location.  The  origin  time  is  said  to 
be  consistent  if  the  difference  between  the  origin  time  of  a candidate  detection 
bulletin  and  that  of  the  top  ranked  'key'  detection  bulletin  is  less  than  the  ex- 
pected timing  error  implied  by  the  distance  between  the  grid  points  of  potential 
event  locations.  To  improve  the  precision  of  the  location  search,  the  single- 
site module  association  tests  are  repeated,  in  toto,  several  times.  At  each 
stage,  the  output  location  estimate  is  taken  as  a new  initial  location.  The  grid 
size  is  halved  and  the  search  is  repeated  to  find  the  maximum  number  of  assoc- 
iates with  consistent  origin  times.  This  is  repeated  until  a prescribed  precis- 
ion is  attained  for  the  depth  constrained  location  estimate.  After  completing 
the  single-site  association  tests,  a depth  search  is  performed  to  find  the  focal 
depth  with  minimum  variance  error  in  predicted  station  arrival  times.  This 
search  is  also  performed  successively  on  a grid  of  depths.  The  search  is 
carried  out  to  the  desired  level  of  precision  by  successive  binary  partitions  of 
the  grid  distance  between  the  preceding  best  depth  estimate  and  neighboring 
grid  points.  Finally,  all  of  the  arrival  times  of  associated  detection  bulletins 
are  used  for  a linearized  least  squares  location  estimate  using  the  method  of 
Geiger  (19101,  further  developed  by  Bolt  (I960)  and  Flinn  (1965). 

A mathematical  description  of  the  Kalman  Filter  and  Geiger's 
method  is  given  in  Appendix  A. 
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The  rationale  for  using  the  Kalman  Filter  or  collapsing  grid 
method  to  perform  the  DAP  is  to  develop  a general  systematic  procedure  for 
selecting  detections  associated  with  an  event  out  of  a detection  bulletin  list 
consisting  mostly  of  false  alarms  from  noise,  coda  of  preceding  events,  and 
later  phases  of  preceding  events.  The  main  difficulty  is  to  be  able  to  assoc- 
iate and  locate  smaller  events  shadowed  by  much  larger  events  occurring  about 
the  same  time.  The  Kalman  Filter  provides  a mathematical  basis  for  obtaining 
sequential  estimates  of  location  and  location  errors.  These  can  be  used  as  a 
basis  of  association  testing.  The  collapsing  grid  can  be  used  to  efficiently 
weed  out  false  alarms.  These  inadvertently  get  into  the  associated  set  of  de- 
tections causing  false  event  declearations  and  also  causing  the  possible  loss 
of  bulletins  which  otherwise  would  have  been  correctly  associated.  The  al- 
gorithms used  to  perform  the  previously  described  sub-system  functions  still 
need  to  be  improved  to  operate  satisfactorily  at  high  false  alarm  rates  mostly 
caused  by  event  interference. 
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SECTION  IV 

EVALUATION  AND  RESULTS 

The  DAP  was  tested  on  four  flays  ol  simulated  network  data. 
This  test  used  the  network  simulator,  which  has  been  reported  upon  earlier 
(Shoup  and  Sax,  1974).  The  network  simulator  is  a model  of  earth  seismicity 
attached  to  a network  of  simulated  automatic  detectors.  The  time  period  in- 
cluded 100  simulated  events.  The  array,  single-site,  and  mixed  modes  were 
tested.  Since  the  network  simulator  assumes  an  array  network,  the  single- 
site  simulation  required  that  the  array  detection  bulletins  be  reduced  to  equi- 
valent single-site  detection  bulletins.  This  was  achieved  by  having  the  DAP 
ignore  the  array  derived  data,  i.e.  , azimuth  and  dT/dA,  for  those  stations 
specified  as  single-  sit  os. 

The  simulated  detectors  generated  false  detection  bulletins  at 
a rate  of  12  per  day.  Interference  problems  such  as  large  timing  errors,  and 
detection  errors  caused  by  signal  coda  and  later  phases  were  also  simulated. 

The  reason  for  using  a false  alarm  rate  from  ambient  noise  of 
12  per  day  was  that  Shoup  and  Sax  (1974)  in  a previous  report  on  the  DAP  found 
that  this  was  an  optimum  false  alarm  rate  for  operating  automatic  detectors. 
Higher  false  alarm  rates  of  detectors  degraded  performance  by  tending  to  sat- 
urate the  bu  letin  buffer.  Lower  false  alarm  rates  made  it  more  difficult  to 
obtain  the  four  bulletins  required  to  locate  smaller  events.  Obviously,  as  the 
DAP  algorithms  are  improved,  then  operation  at  higher  false  alarm  rates  be- 
comes feasible.  The  interaction  between  station  detector  threshold  strategy 
and  DAP  performance  is  a very  important  factor  in  the  overall  detection  per- 
formance of  world-wide  seismic  networks  and  needs  more  careful  investigat- 
ion. 
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Before  examining  results,  it  is  useful  to  look  at  where  things 
can  go  wrong.  Numerous  bulletins  originating  from,  but  not  correctly  asso- 
ciated with  a large  event  can  prevent  overshadowed  smaller  events  from  .at- 
taining top  ranked  status  in  the  bulletin  buffer.  In  the  configuration  shown  in 
Figure  III  - 1 , events  occurring  simultaneously  must  be  processed  serially. 

Chit  of  such  multiple  occurrences,  the  smaller  events  can  be  missed;  even 
with  a very  large  number  of  bulletins  reporting. 

For  example,  if  several  events  go  off  at  about  the  same  time, 
the  largest  event  detections  will  be  at  the  top  of  the  /.-ranked  workspace  of  the 
bulletin  buffer.  After  that  event  is  associated  and  located,  all  bulletins  with 
small  arrival  time  errors  are  pulled  out  of  the  workspace.  Possibly  a substan- 
tial number  of  later  phases  and  coda  detections  w'ill  still  be  left  in  the  work- 
space. Some  of  these  will  graduate  to  the  top  of  the  workspace  and  atte  mpts 
will  be  made  to  associate  them  with  other  detections  in  the  workspace.  Obvi- 
ously since  they  are  not  P-waves  of  an  event,  they  arc  unlikely  to  associate. 
They  tend  to  clog  the  bulletin  buffer  during  which  period  of  time,  good  P-wave 
detection  bulletins  will  be  aged  out.  One  way  to  avoid  this  problem  is  to  pro- 
vide the  capability  of  performing  association  tests  on  many  key  detections  in 
parallel.  If  a detection  ages  out  on  trying  to  associate  with  a key  detection  at 
one  level,  it  is  pushed  down  a level  and  tries  to  associate  with  a key  detection 
at  the  next  level.  This  could  possibly  prevent  large  false  detections  from  clog- 
ging the  workspace.  The  DAP  strategies  considered  in  this  report  do  not  ade- 
quately handle  this  problem  and  therefore,  a substantial  number  of  small  events 
following  larger  events  failed  to  be  located;  in  fact,  one  event  with  15  good  de- 
tections was  missed  possibly  for  this  reason. 

Decision  errors  in  association  tests  can  cause  large  errors  in 
event  locations  obtained  by  successive  use  of  the  Kalman  Filter.  The  operation 
of  the  Kalman  Filter  is  a bootstrap  operation.  In  association  testing,  since  the 
probability  of  a coincidental  overlapping  of  an  event  location  by  the  confidence 
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region  of  a small  array  location  error  ellipse  is  on  the  order  of  0.  01,  there 
is  a small  but  yet  significant  probability  for  a false  alarm  to  be  accepted  into 
a set  of  associated  detections.  Since  the  false  alarm  is  likely  to  have  a large 
arrival  time  error,  any  subsequent  location  is  more  likely  to  exceed  the  two 
degree  requi rement  for  location  accuracy  than  would  otherwise  be  the  case. 
Presently,  there  is  no  means  of  extracting  these  erroneously  associated  bul- 
letins. This  can  cause  missed  associations  and  substantial  errors  in  the  final 
location.  The  single- site  association  method  using  exhaustive  search  can  re- 
cover from  such  association  errors.  Erroneously  associated  bulletins  can  be 
eliminated  as  the  grid  length  is  successively  halved.  The  single- site  associa- 
tion test  may  occasionally  miss  bulletins  by  incomplete  coverage  of  all  of  the 
space  required  to  capture  all  possible  event  locations.  If  necessary,  this 
could  be  prevented  through  overlapping  the  event  coverage  of  single- site  as- 
sociation tests,  by  increasing  the  threshold  of  the  origin  time  consistency  test. 

Table  IV- 1 shows  the  results  of  locating  the  events.  The  most 
significant  information  is  given  by  the  success  or  failure  to  locate,  given  that 
four  or  more  accurate  detection  bulletins  reported.  The  best  results  were 
given  by  the  single-site  associator.  Successes  are  defined  as  location  errors 
less  than  200  kilometers  and  origin  time  errors  less  than  15  seconds  for  those 
events  of  which  four  or  more  stations  detected  the  arrival  time  of  P-waves. 
Failures  are  defined  as  errors  exceeding  those  limits  or  failures  to  associate 
for  those  events  of  which  four  or  more  stations  detected  the  arrival  time  of 
P-wave  signals.  Thus,  worse  results  were  obtained  with  an  array  network 
despite  the  availability  of  the  additional  array  measurement  information  on 
array  detection  bulletins.  This  clearly  indicates  serious  problems  with  the 
association  test  of  the  array  module.  The  importance  of  not  missing  proper 
associations  is  demonstrated  by  the  mixed  network  results.  The  association 
tests  rejection  threshold  were  set  for  error  ellipses  greater  than  la,  2a, 
and  3a  (the  single- site  and  array- site  networks  used  3a  association  tests, 
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TABLE  IV-1 


DETECTION  ASSOCIATION  PROCESSOR  RESULTS 


Single-Site 

A r ray 

Mixed  Network 

1(7 

Rejection 

2(7 

Rejection 

3a 

Rejection 

Total  No. 

of  Events 

128 

100 

100 

102 

102 

Successes 

40 

i 

30 

12 

26 

29 

Failures 

17 

17 

35 

24 

21 

Interfe  red 

Events* 

29 

21 

21 

22 

22 

Events  De- 
tected At 
Less  Than 

4 Stations 

42 

32 

32 

30 

30 

Interfered  events  are  failures  of  the  automatic  detector  to  properly 
time  events.  Large  timing  errors  were  caused  by  misidentification 
of  the  P phase  due  to  interference  with  other  event  phases.  Most  of 
these  large  timing  errors  could  easily  be  prevented  by  improving 
the  logic  of  the  simulated  automatic  detector  used  for  this  analysis. 
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whe-re-  CT  is  de'finetl  as  the  mean  standard  deviation  of  tin-  error  ellipse).  Note 
the  eatestrophie  drop  in  the  success  ratio  for  the-  2o  and  lo  re  lection  threshedds. 
This  is  due  to  not  discarding  larger  detections  associated  with  located  events. 

The  success  ratio  for  the  mixed  network  was  less  than  array- 
site  network  and  the  single'- site  network.  This  suggests  the  possibility  that 
the  interface  between  the  array  module  and  single-site  module  may  not  be  pro- 
perly configured. 

Table  IV- 2 shows  that  the  detect  ion  association  processors 
location  errors  appear  to  be  consistent  with  correct  assoc  iation  and  hypocen- 
ter  location.  An  interesting  observation  is  that  the  origin  time'  error  and  depth 
error  standard  deviation  were  halved  for  the  array  module.  Looking  back  at 
Figure  111-1,  note  that  the-  depth  grid  search  was  performed  by  the  array  mod- 
ule- but  not  by  the  single-site*  module'.  The'  dotted  line-  indicates  that  the-  depth 
search  was  not  ye't  implemented  into  the  single'- site'  module.  Thus  it  appears 
that  a depth  grid  search  for  minimum  arrival  time  errors  can  lead  to  large  re- 
duction in  depth  and  origin  time-  estimation  error.  Since-  we  might  otlu-rwise 
have  not  se'e-n  this  re'sult,  it  was  perhaps  fortunate-  that  tlu-  de-pth  search  was 
not  put  into  the-  single'- site  module. 

Figure-  IV- l shows  the  re-sults  obtaine'd  for  the-  c'vents  analyzed 
by  the  single-site'  module'.  Note'  the-  large'  number  of  failures  with  four  or 
more  bullt'tins  at  magnitudes  from  4.3  to  4.').  This  suggests,  strongly,  ser- 
ious proble-ms  with  re-al-time-  processor  servicing  the-  bulle'tin  buffer.  Some- 
possible  source's  of  the-se  problems  we're'  prc*viously  me-ntione-d.  By  going 
back  to  the'  bulletin  list  and  checking  these  failures,  most  were  found  to  be- 
pre-ce'e-de'd  by  large'!'  e-vents.  The'  seriousness  of  this  proble'm  is  further 
shown  by  Table'  IV-3.  Note  that  e've'nts  were  missed  with  as  many  as  l1'  valid 
de'te'ction  hulle'tins.  Assuming  that  the  4 and  k bulle'tin  cases  ol  failing  to  lee- 
cate-  the*  e've'nt  are-  as  expe'ete'd;  those'  e've'nts  with  faile'd  location  e-stimate's, 
given  (>  or  more-  hulle'tins,  appear  to  be-  anomalous,  t'n  that  basis,  it  is  e-sti- 
mated  that  75"i  of  the  errors  could  be'  e'liminate'd  by  modifying  tlu-  re'al-time- 
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TABLE  IV- 2 


LOCATION  ERROR  ANALYSIS 


95%  Error 
Ellipse 
(km) 


1 2 2 

-r.L(e.  + e 

N lat  Ion 


Single-Site 

A r ray 

Mixed  Network 

1<T 

Rejection 

2 a 

Rejection 

3 a 

Rejection 

152.  06 

214.  7 

224. 01 

1 39. 32 

1 39.  06 

1. 08 

-0.  54 

0.  16 

1.  72 

-0.  94 

6.  61 

3.  03 

20.. 85 

9.  61 

4.  84 

13.  86 

-4.  69 

35.  96 

12.  55 

-0.  87 

60.  73 

36.  57 

158. 63 

77.  37 

46.  28 

21.  35 

25.  4 

25.  9 

20.  4 

20.  4 

variance  of  location  errors  in  kilometers 

mean  timing  error 

standard  deviation  or  origin  time 

mean  depth  error 

standard  deviation  of  depth  error 
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TABLE  IV -3 

ANALYSIS  OF  LOCATION  FAILURES 


Nun.Ler  of  Detected 

P- Wave  Sigr'Us 

Magnitudes  ot 

Missed  Events 

4 

4.  0,  4.  0,  4.  2,  4.  3,  4.  4 

5 

4.  3 

6 

4.  4,  4.  5 

7 

4.  3,  4.  5 

8 

4.  3 

9 

4.  7 

10 

4.  4 

11 

4.  4 

12 

4.  5 

13 

4.  9 

14 

15 

4.  3 

'****■*'  •■+>*&** 


prov'mtior  wbu  b m-rvii  p*  thu  ballot  in  buffor.  It  is  ooniooturod  that  suoh  .1 
niodlu.it  ion  ooubl  rottult  in  .1  sun  oss  to  failuro  ratio  of  about  It'  to  I. 
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The  detection  association  processor  described  in  this  report 
indicated  a capability  of  associating  and  locating  about  70'/'  of  the  events  with 
four  or  more  satisfactory  P-wave  detections.  This  was  done  on  a completely 
automated  basis.  It  can  handle  the  detection  association  and  location  function 
for  single  sensor  networks,  array  networks,  or  mixed  networks  containing 
both  types  of  stations.  The  accuracy  of  the  locations  appeared  to  be  satisfac- 
tory and  comparable  with  the  visual  linearized  and  iterative  least  squares 
hypocenter  determination.  In  fact,  a depth  gird  search  for  least  squares  ar- 
rival times  produced  depth  determinations  which  were  almost  a factor  of  two 
more  accurate  than  those  hypocenter  determinations  that  were  unconstrained 
with  respect  to  depth  of  focus. 

The  results  clearly  indicated  that  more  work  is  needed  to  im- 
prove the  detection  association  processor.  The  main  discrepancy  is  the  lack 
of  a satisfactory  criteria  for  removing  aged  and  redundant  detections  from 
the  workspace  of  the  processor.  For  example,  later  phases  should  be  identi- 
fied and  removed  from  the  detection  association  processor  workspace  as  soon 
as  possible  after  an  event  location  is  declared.  Also,  the  process  by  which 
detections  are  aged  and  discarded  must  guarentee  that  the  detection  associa- 
tion tests  be  performed  on  a sufficient  number  of  key  detections  to  associate 
and  separate  several  events  occurring  simultaneously  over  a reasonably  wide 
range  of  magnitude. 

The  association  process  needs  to  be  improved.  The  initializa- 
tion of  key  detections  which  are  association  tested  against  other  detections 
needs  to  be  handled  as  a parallel  process  rather  than  a serial  process.  In 
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other  words,  several  events  should  be  associated  simultaneously  rather  than 
one  at  a time.  Also,  the  tests  used  to  perform  the  association  of  detections 
need  to  he  improved.  Time  difference  constraints  should  be  imposed  on  all 
possible  pairs  of  association  candidates  based  on  the  time-distance  relation- 
ship for  P-wavcs.  Other  criteria  might  be  considered  in  selecting  key  detec- 
tions other  than  the  magnitude  of  the  detector  output.  For  example,  one  such 
criteria  is  the  time  elapsed  at  a station  from  the  previous  larger  magnitude 
detection;  another,  arrival  times  earlier  than  the  current  key  detection.  The 
first  time  criteria  assures  a lower  probability  of  missed  associations  due  to 
preceding  larger  events;  the  second,  that  older  detections  be  tried  as  key 
detections  before  being  timed  out. 

The  suggested  modifications  of  the  automatic  detection  associ- 
ation processor  should  build  upon  other  association  processors  and  location 
procedures  now  in  use,  as  described  by  Der  (1977),  Flinn  ( 1 9t> 5 ) , and  Rolt 
(I960).  Although  the  performance  shown  here  is  less  than  that  of  existing 
association  processors  which  are  80%  to  90%  efficient  in  associating  and  lo- 
cating events  compared  to  the  70%  efficiency  demonstrated  here,  it  is  probable 
that  performance  much  closer  to  100%  efficiency  might  be  possible  after  work- 
ing out  obvious  faults  in  the  program. 
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APPENDIX  A 


A mathematical  description  of  the  Kalman  Filter  and  Geiger's 
method  is  given  to  specify  the  detailed  methodology  of  the  DAP  algorithms 
used  in  this  report. 


Kalman  Filter 

Suppose  there  is  a non-linear  system 


where 


Define 


and 


/.  = h(x)  + V 


(A  - 1 ) 


7.  - (o  x 1)  measurement  variable  matrix, 

x = (n  x I)  state  variable  matrix, 

h(x)  = (p  x 1)  non-linear  function  matrix,  and 

V = (p  x 1)  measurement  error  matrix. 


E(VVT| 


R,  a (p  x pi  positive  matrix. 


(A- 2) 


E (x  - x)  (x  - x)T  j = 


M. 


(n  x n)  matrix,  error  co- 
variance  before  measurement 


A - 1 


(A-  3> 


where 


x = the  predicted  state  variables  before  measurement. 

We  can  find  a controlling  function  J(x)  which  is  often  called  the  performance 
criteria,  or  the  performance  index. 

J(x)  = ij(x-3f)TM  1 (x  - x)  + |z  - h(x)jT  R ]^z-h(x)j|.  (A-4) 


Minimizing  J(x)  with  respect  to  x obtains  an  estimate  of  x 


x = x + K(x) 


- h(x)J 


where 


K(x) 

C(x) 

P 


T - 1 

PG  R , a (n  x p)  Kalman  Filter  matrix, 
3h(x)| 

a"  - , (pxn)  matrix, 

ffx  I x=x 

E £(x  - x)  (x  - x)T]  , (n  x n)  matrix,  error 
covariance  after  measurement, 

.t1 


(A-  5) 

(A-  6) 
(A  — 7) 


t r t! - i 

M - MG  R + GMG  J 


GM. 


(A-  8) 


In  our  system,  we  have  three  state  variables  and  three  mea- 
surement variables.  Their  relationships  are  given  below; 


Z, 


Z„ 


h j (x) 
h2(x) 
h3(x) 


= cos  (A) 

= sin(A)  sin(Az) 

= TA 

= sin(xE)  sin(Xg)  + cos(xE)  cos(xg)  cos(yE-yg) 
= cos(xE)  sin(yE-ys),  and 


T + T 
o 


(A -9) 


( A - 1 0 ) 


A-  2 


r 


! 


who  re 


= epicentral  distance  between  the  event  and  the  station, 
- the  azimuth  angle  of  the  observer, 

= the  latitude  of  the  event, 

= the  longitude  of  the  event, 

= the  latitude  of  the  station, 

= the  longitude  of  the  station, 

= the  origin  time  of  the  event, 

= the  travel  time  of  the  event,  and 
= the  arrival  time  of  the  event. 


The  R Matrix  : 


R is  a (p  x p)  square  matrix.  It  is  obtained  from  measure- 
ments, and  is  usually  a diagonal  matrix  independent  of  state  variables.  It 
plays  an  important  role  in  iterative  computations  of  the  event  location. 


Differentiating  equation  (A -9)  we  have 


V = dZ  = - sin(A ) ~r~  dp,  and 
i 1 op 


(A  - 11) 


where 


V = d7.  = cos(  A)sin(Az)  — — dp  + sin(A  )cos(Az)d(Az),  (A -12) 

Z Z op 


dT 

dA  ‘ 


The  error  of  V includes  a measurement  error  e(T  ) which  was  set  equal 

J t\ 

to  4 sec,  and  an  error  due  to  imperfect  location.  Therefore, 


v = 4 + -f—  dp 

3 9A  3p  p 


(A-13) 


The  change  in  p in  equations  (A - 11)  and  (A -12)  is  identical  and  implicit  from 
the  definition  of  and  7,^.  13ut  since  the  arrival  time  errors,  dZ^,  are 


uncorrelated  with  dZj  and  dZ^,  the  change  in  p is  uncorrelated  in  equations 
( A - 11 ) , ( A - 1 3 ).  and  equations  (A -12)  and  (A -13). 

The  expected  operator  E is  a linear  operator.  If 
R - E £(x  - x)  (x  - x)^]  , it  operates  on  each  term  of  the  matrix.  Therefore, 


R11 

= EIVjVj) 

= 

[ s in ( A ) dp]2  , 

R12 

= EtVjV^ 

= 

-£  sin(2A)  sin(Az)  dp]  2 , 

R13 

= E(VjV3) 

= 

0 , 

R21 

= E(V2Vl) 

= 

R12  ’ 

R22 

= E(V2V2) 

= 

^ cos(A)  sin(Az)  dp]  2 , 

+ 

£ sin(A)  cos(Az)  d(Az)]  2 , 

R23 

= E<V2V3) 

= 

o , 

R31 

= E(V3Vj) 

= 

0 , 

R32 

= E(V3V2> 

= 

0 , and 

r,„ 

33 

= E<V3V3> 

= 

r dT  d±  , 12 

Op  dpJ  • 

The  G Matrix: 

C(x) 

is  a (p  x n)  matrix. 

Each  element  can  be  derived  from 

equation  (A-10),  so, 

C11 

tfh^x) 

**E 

= 

cos(xF)  sin(xs) 

- 

sin(xE)  cos(xs)  cos(yE-ys)  , 

Shj  (x) 

°2i  2 ~a^ 

9h^(x) 

C22  = ~~^E 

9h3(x) 
C23  = ~T^ 


-cos(xE)  cos(xg)  sin(yE*ys)  , 

cos(xE)  cos(yE-yg) 

5h3(x)  ah^x)  /ahj(x) 
dxE  dVE  / dxE 


A - 5 


X 


I 


the  state  variables,  representing  the  latitude  (x  ), 

El 

longitude  (y^.),  origin  time  (T^),  and  depth  (z)  of 
an  event  , and 


Jk 


who  re 


h(x)  = T + T(A, z)  . 
o 


(A  -1 5) 


A = the  epicent  ral  distance  between  the  event  and  the 

observer,  and 

V = the  measurement  error. 

Linearizing  h(x)  about  an  initial  estimate  x - x we  have 


Ta  = h^>  + | x _ (X  . X)  + 0 (|x- xT)  + V . (A -16) 

For  small  | x - x | , 


ta  ’ 

h(x)+ 

$h(x)  , .. 

9 * (x  - x)  * V , 

(A -17) 

R 

CA 

(A-18) 

R 

ta- 

h(x>  - V , 

C 

0h(x 
8 x 

^ 1 — , and 

1 x=x 

(A -19) 

Suppose  there  are  n stations  or  n measurements,  then  R is 
an  (n  x 1 ) row  vector;  C,  an  (nx4)  matrix;  and  A,  a (4x1)  row  vector.  At  each 
station,  C has  four  elements  and  their  expressions  are; 


A-b 


i: 

cu  - ( 

'«h(x)\ 

= 

(e\ 

0A  \ 

. i 1.2... 

i 

cu  1 

f 

= 

(V. 

\d± 

*yE)i 

. i 1,2... 

i 

Ci3  - 1 

s 

l 

. i 1.2... 

i 

Ci4  - | 

f 0h(x)\ 

^ 9*)i 

s 

(*) 

i 

. i 1,2... 

• . n ; 


(A-  20) 
, n ; a ml 


Minimizing  the  norm  of  error 

II  II 2 


we  have 


1.  r.  , 


(R  - CA)  (R  - CA) 


T - 1 T 
(C  C)  C R 


x = x + (CTC)'1  C1  R . 


(A- 21) 


This  calculation  can  be  continued  iteratively  until  tin-  term 


T - 1 T 

(C  C)  C R becomes  insignificant. 


i. 


